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ABSTRACT 



A systematic, rigorous, and complete investigation of the Bloch equations in time-harmonic 
driving classical field is performed. Our treatment is unique in that it takes full advantage 
of the partial fraction decomposition over real number field, which makes it possible to 
find and classify all analytic solutions. Torrey's analytic solution in the form of exponen- 
tially damped harmonic oscillations [Phys. Rev. 76, 1059 (1949)] is found to dominate 
the parameter space, which justifies its use at numerous occasions in magnetic resonance 
and in quantum optics of atoms, molecules, and quantum dots. The unorthodox solutions 
of the Bloch equations, which do not have the form of exponentially damped harmonic 
oscillations, are confined to rather small detunings 5 2 < (7 — 7t) 2 /27 and small field 
strengths f2 2 < 8(7 — ■y t ) 2 /27, where 7 and 74 describe decay rates of the excited state 
(the total population relaxation rate) and of the coherence, respectively. The unorthodox 
solutions being readily accessible experimentally are characterized by rather featureless 
time dependence. 
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1 Introduction 



The Bloch equations in time-harmonic driving classical field have been employed over 
several decades as an important tool in studies of many different physical phenomena 
Pai^llSllSISllSllZllSllSl Hia IIH H3l H3l HSl [TO]- The equations [see © below] underline 
the theory of magnetic resonance [U 121 El El E] and the quantum optics of a two-level atom 
(molecule, spin, ion, etc) driven by a classical field [31 El El El [91 [101 [Til H21 [13l [H] . The 
latter problem is one of the most discussed and is at the heart of the theory of self-induced 
transparency [HI E], the susceptibility of an ensemble of atoms (spins, ions, etc) in gain 
media [9], and a number of other optical phenomena [SI El HOI [HI EE21 EES]- Recently the 
problem has been extensively studied in connection with the proposed use of atoms and 
molecules as a triggered single-photon emitter [Til 02], a single-photon emission of reso- 
nantly driven semiconductor quantum dot in a microcavity [15J, and decoherence in dc 
SQUID phase qubits [IB] . Advances in the fabrication of single defect centres in diamond 
enable one to investigate two-level diamond-based single-photon emitters at room tem- 
peratures [T71 [18] . At the same time the single defect centres enable nanoscale magnetic 
sensing, and hence a nanoscale imaging magnetometry, with an individual electronic spin 
in diamond under ambient conditions [191 EH [2Tj - 

The Bloch equations form a linear system of three ordinary differential equations [see 
(J2J) below], which can be formally solved in the following three steps: 

• (i) applying the Laplace transform, whereby the system of differential equations 
reduces to a linear algebra problem; 

• (ii) solving the ensuing linear algebra problem (for instance, by means of Cramer's 
rule) ; 

• (iii) applying the inverse Laplace transform. 

That was also the original Torrey's approach [21 E], who established a general time 
dependence of each Bloch variable ui, I = 1, 2, 3, in the form of exponentially damped 
harmonic oscillations 

ui(t) = A + A ie " Kl * + A 2 e~ bt cos(st) + {A^/s)e~ bt sm(st). (1) 

In order to determine explicit analytic solutions, essential to Torrey's approach was to find 
out a negative real root —K\ (it always exists - see section E]) of a characteristic determinant 
A [see (fTDT) below] of the Bloch equations. Assuming the knowledge of Ki, Torrey [2] 
determined the constants Aq, A\ as certain limits of Laplace transforms [see (|55|) . (1351) 
below], and determined the remaining constants A 2 , A3 using the initial conditions and 
the knowledge of A and A\ (see Eqs. (45)-(47) of [2]). The root parameters b and s were 
then determined on comparing expansion coefficients of the cubic polynomial A(p) in the 
Laplace transform variable p against corresponding expressions of the coefficients in terms 
of cubic roots [see Eq. (48) of [2] and ffTB"]) below]. Not aware of Cardano's formula [21] - a 
jewel of the mathematics of the 16th century [25] - Torrey managed to formulate explicit 
analytic solutions merely in three special situations [2]: 
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• strong collisions: 7 = 7t (Tl); 

• exact resonance: 5 = (T2); 

• intense external field: fl/jt ^ 1 (T3). 

Additionally, it had long remained unnoticed that Torrey's general solution ([T]) had been 
confined to the parameter range when the discriminant D [see (fT8l) below] [27] of the 
characteristic determinant A is negative (see section EJ). Indeed, general solutions in the 
parameter range D > are nonos dilating and can be represented as a sum of three 
exponentially damped terms, wherein each term corresponds to a real cubic root of the 
characteristic determinant A (see section [S]). 

Surprisingly enough, the limitation of Torrey's general solution ([1]) to the parameter 
range D < and the unawareness of Cardano 's formula had been since without exception 
repeated in the literature [21 [7J and in textbooks (cf section 3.5 of Ref. [Sj) for several 
decades. The latter is quite remarkable given the role that Torrey's analytic solutions 
played in quantum optics and magnetic resonance literature [7J [8] . It took several decades 
before Hore and McLauchlan had finally made use of Cardano's formula [21] and formally 
determined the cubic roots [I]. However, their analytic solution has other deficiencies to 
be discussed in section [SJ 

In what follows, we closely follow the original Torrey's approach [21 HI E] and supple- 
ment it with two additional elementary tools: 

• (i) the partial fraction expansion, or partial fraction decomposition (PFD) [23J over 
the field of real numbers (see Appendix |A]) , combined with the classification of roots 
of a real cubic polynomial, and 

• (ii) a fully compensated formula for real roots of a real cubic polynomial [see (|38|) 
below], which has recently been derived [2S] from Cardano's formula |24j . 

The first tool will enable one to classify and, up to an explicit knowledge of the roots of 
A, determine all the possible functional forms of solutions of the Bloch equations ([2]) for 
all values of D = 0. The second tool will allow one to determine the roots of A, and hence 
solutions of the Bloch equations, explicitly. 

The outline of the article is as follows. In section [2] we summarize notation and give 
some basic definitions. In section [3] a full advantage of the partial fraction decomposition 
over real number field is made, all possible solution types (e.g. of non Torrey type) are 
classified [see ( JTT)) ]. In section @] the parameter range of unorthodox solutions is deter- 
mined. In section [5] decay constants corresponding to each solution type are determined 
by Cardano's formula. In section El sufficient and necessary conditions for doubly and 
triply degenerate real roots are established in terms of the model parameters. Steady 
state solutions and remaining numerical constants of different solution types [see flT7|) ] are 
determined in section [7J In section [8] connection to earlier results and various conditions 
of the applicability of our results are discussed. We then conclude by section [91 A number 
of formulae and intermediary steps are relegated to Appendices \M and [Bl 
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2 Bloch equations 



The Bloch equations are a linear system of differential equations with constant coefficients 
for a Bloch vector with components (u, v, w), also known as the Bloch variables [H El IH 

El Eld Eli, 

v! = —jt u — 8v, 

v' = — jtv + 5u + flw, 

w' = —jw — Qv + •yweq, (2) 

where the prime denotes derivative with respect to a parameter t defined in Table 1, 
and w eq is an equilibrium value of w. Occasionally we denote the Bloch variables as 

(u,v,w) = (ui,u 2 ,u 3 ). 

Table 1. The Bloch equations parameters in the case of a two- level system and in the 
case of magnetic resonance (MR). 





two-level system 


MR 


t 


time r 


rescaled time gH\T 


7 


l/T x 


V(^iTi) 


It 


1/T 2 


V(gHiT 2 ) 




(e 2 - e x )/K 


Larmor frequency gH 


S 


Uq — U 


(uj - u)/{gH x ) 


Q 


h 


-1 



The meaning of the parameters depends on the problem solved. In the case of a two-level 
system with energies e\ and e 2 > t\ interacting with time-harmonic perturbing classical 
field E{t) = E cos ut with the driving frequency oj and amplitude E , the Bloch variables 
u, v, w are related to the elements of a 2 x 2 hermitian density matrix py = p* t [HE] 

w = 2Repi2, w = 2Impi2, w = p 73 — p 11 . (3) 

The parameter 5 stands for the driving field detuning from the intrinsic resonance fre- 
quency oj , whereas Q, which reduces to the Rabi frequency at zero detuning, accounts 
for the interaction strength - it is entirely determined by the coupling with the driving 
field (e.g. between the electric field amplitude and the transition dipole moment fi) (see 
Table 1). Phenomenological constants 7 and j t describe decay rates of the excited state 
(the total population relaxation rate) and of the coherence, respectively, where T\ is the 
spontaneous emission lifetime and T 2 (typically > J\) is the total dephasing time [HI E]- 
According to (j3J), w is the single atom population difference, also called inversion [HIE]- 
The condition w(0) = —1, which (together with u = v = 0) is often taken as the initial 
condition [H El EH] , means that the two-level atom is initially in its ground state. The 
third equation in (T5]) shows that v is the component effective in coupling to the field to 
produce energy changes. Thus v determines the absorptive (in- quadrature with the field 
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E) component and u the dispersive (in-phase) component of the atomic transition dipole 
moment. Note in passing that w eq is not necessary a thermal equilibrium value, since 
some pump mechanism may be present that causes w eq at equilibrium to have some fixed 
value that is different from its thermal equilibrium value [9]. 

In the theory of magnetic resonance of precessing (nuclear or electronic) spins [TJ |2j HJ 
E], the time constants T2 and T\ correspond to the longitudinal and transverse relaxation 
constants [H El HI E] and, in contrast to the optics, T 2 is always less than or equal to 
Ti [22]. Assuming a constant magnetic field H applied along the z-axis and a time- 
harmonic component 2Hi cos ut applied in the x-direction, the components u, v, w in Eqs. 
([2]) are formed by the respective x, y, z components of a nuclear or electron (macroscopic) 
magnetization M, and w eq is related to the static magnetization M = XoHc where Xo is 
the static susceptibility [TJ |2]. Eqs. d2J) are then valid with the time r together with the 
homogeneous lifetimes T\ and T2 being rescaled by gHi, where g is the absolute value of 
the gyromagnetic ratio (see Table 1). 

Importantly, in any case is the ratio Q/\j — jt\ proportional to the driving field (either 
E or Hi) and (see Table 1) 

x 

(4) 



l7~7t| 





0J - 


CO 




1 


1 






T-2 



3 Classification of possible solution types 

Following Torrey's approach [21 HJ [5] , upon applying the Laplace transform, 

poo 

Kp) = / f(t)e^dt, (5) 



the Bloch equations ([2]) are transformed into the matrix equation 



P + lt 5 
-5 p + lt -n\\v\ = \v Q \, (6) 





( u\ 




( uo \ 




V 


\-[ 










\ w J 



Q p + 7 
where tilde denotes the Laplace transform of the Bloch variables, 

W = w + (7) 
V 

and Uq, Vq, Wo are the initial values of the Bloch variables. Cramer's rule then yields 

where ui, I = 1,2, 3, stands for u, v, and w, respectively, and 

M(P + 7*)(P + 7) + - v 5{p + 7) - W 5Q, I = 1 
fi(p) = { Uo£(P + 7)+Vo(p + 7t)(P + 7) + lW(p + 7t), 1 = 2 ( 9 ) 

-u sn - v n( P + lt ) + w [(p + it) 2 + s% 1 = 3 
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The determinant A(p) of the coefficient matrix is a real cubic polynomial 



A(p) 



+ a 2 p 2 + aip + a , 



(10) 



where 



a 2 = 



a = 



7 + 27t, 

7 2 + 2 77t + o" 2 + fi 2 , 
77 4 2 + 7^ + 7^ 2 - 



(11) 
(12) 

(13) 



The multiplication of both the numerator and denominator on the right-hand side of 
(jHJ) by p changes each fi(p) defined by Eqs. (jSJ) into a cubic polynomial, whereas the ratio 
on the right-hand side of (jHJ) becomes the ratio of a cubic and quartic polynomials. Given 
that A(p) is a rea/ cubic polynomial [see (TTTi) to (TT3T) for the coefficients a, in (flOl) ]. this 
suggests the application of the PFD over the field of real numbers [see fj67j) in Appendix 
lAj . The PFD enables one to express the fractions such as that in (jHJ) as a sum of much 
simpler fractions, whose inverse Laplace transform may be readily available. To this end 
note [e.g. by the fundamental theorem of algebra (I66p ] that as any real cubic polynomial, 
A(p) can have either (i) one real root and a pair of complex conjugate roots or (ii) three 
real roots. Given that all the coefficients aj in (fTUj) are positive real numbers, one can 
prove additionally that: 

• A(p) has always at least one negative real root (PI); 

• if A(p) has three real roots, they are all negative (P2); 

• if not all roots of A(p) are real than there is one real root and two complex conjugate 
(c.c.) roots (P3); 

• if at least two roots of A(p) coincide, they are all negative real roots (P4). It may be 
that A(p) has a double real root and another distinct single real root; alternatively, 
all three roots of A(p) coincide yielding a triple real root. 

The properties PI and P2 can be shown to be a straightforward consequence of the 
intermediate value theorem when applied to a cubic polynomial A(p) with positive real 
coefficients aj. Obviously, a sufficient condition for PI is A(0) = ao > 0, whereas a 
sufficient condition for P2 is 



The property P3 follows from PI and the fundamental theorem of algebra (I66p applied to 
a real cubic polynomial. Eventually, the property P4 follows upon combining PI to P3. 
Note in passing that we disregarded a special case of ao = 0, which is treated at the end of 
section To this case belongs also the trivial undamped case 7 = 7 t = (which implies 
both ao = and a 2 = 0), in which case the Bloch equations ([2]) describe the precession of 
a classical gyromagnetic moment in a magnetic field [21 IB] ■ 



A(p) > 



(p>0). 



(14) 
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Given the properties PI to P4, and on writing negative real roots of A(p) 
j = 1,2,3 (kj > 0), the application of the PFD [see (157]) in Appendix IA~] enables one to 
decompose each /j(p)/A(p) as 



J(p) = P/(P) 
A(p) pA(p) 



_|_ 

X _|_ 
X _|_ 

X I 
p 



Ai 

Al 

M 



A 2 {p+b)+A 3 
(p+b) 2 +s 2 ' 
A 2 _,_ A 3 

A 2 



P+K2 
M 



+ 



P+K3 ' 

As 



P+Kl (P+Kl) 2 



hK 2 ) 2 ' 

Aa 

O+Kl) 3 ' 



(pair of c.c. roots) 
(distinct real roots) 
(double real root) 
(triple real root). 



(15) 



Herein and below the index I labeling different Bloch variables will be suppressed unless 
explicitly required. In the first line of Eqs. ( 1T5]) . corresponding to cubic roots — K\, z, and 



z, we have recast the irreducible quadratic factor (p— z)(p— z) in Torrey's form (p+b) 2 + s 2 . 
Given that the coefficients a 2 and a in Eqs. ( ITT]) and flTB"]) can be alternatively expressed 
in terms of cubic roots Zi, z 2 , and z 3 of A(p) as a 2 = — (zi + z 2 + z 3 ) and a = —ziz 2 z 3 , 
the Torrey constants b and s can be entirely expressed in terms of K\, a 2 and a , 



-Re z 



a 2 - Ki 



VP 



& 2 



(16) 



On substituting ( TT5l) back into (JHJ) one can perform the inverse Laplace transform (see 
Appendix [B]). Note in passing that the property P2 additionally guarantees the very 
existence of the inverse Laplace transform. (Indeed, if one of the real roots were positive, 
one would face singular integrals.) After the inverse Laplace transform one arrives at the 
following complete set of the possible functional forms of solutions of the Bloch equations 



ui(t) 



A + A x e~ nYt + A 2 e' bt cos(st) + (A 3 /s)e' 
A + A x e~ Kxt + A 2 e~ K2t + A 3 e~ Kst , 
A + A x e~ Kxt + A 2 e~ K2t + A 3 te~ K ' 2t , 



bt 



A + A x e 



-Ki t 



A 2 te~ Klt + ±A 3 t 2 e~ 



Kl t 



sin(st), (pair of c.c. roots) 
(distinct real roots) 
(double root) 
(triple root). 



(17) 

The PFD's given by Eqs. ( 1TB1) exhaust all the possible PFD's of the ratio f(p) / A(p). There 
is no other decomposition possible. Only the first line in Eqs. ( fTBT) and (ITT]) corresponds 
to the Torrey solution [2], whereas remaining lines yield unorthodox (i.e. non Torrey) 
solutions. Importantly, the PFD guarantees that the respective sets of numerical constants 
Aq, Ai, A 2 , and A 3 for different Bloch variables, together with constants b and s, are all 
real numbers. 



4 Parameter range of unorthodox solutions 

To this end, we have not verified yet if any of the above solution types could be attained 
within the physical range of parameters in the Bloch equations ([2]). In order to proceed, 
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one needs to determine the value of the discriminant D of A(p). According to 

L> = a 1 a 2 — 4a^ — 4a a2 — 27a + 18a aia2 

= ( 7 - lt ftf (ft 2 + 205 2 ) - 4( 7 - 7t ) 2 5 2 [(7 - It? + 25 2 ] 

-4(5 2 + ft 2 ) 3 . (18) 

In terms of the roots, the discriminant of a general polynomial of the nth order is given 
by W\ 

D = af - 2 J] (* - ^ (19) 

where a n is the leading coefficient [a n = 03 = 1 in (flUj) ] and Zi,...,z n are the roots 
(counting multiplicity) of the polynomial. For a cubic polynomial with real coefficients, 
the use of the general form of the discriminant [27] together with its real value [see f)18p ] 
enables one to relate the nature of roots to the value of D as follows: 

• D < 0: a cubic polynomial has one real root and two complex conjugate roots; 

• D > 0: there are three distinct real roots; 

• D — 0: at least two roots coincide, and they are all real. 

To this end one can show that Torrey special solutions labeled by Tl and T3 introduced 
in section [1] are always described by exponentially damped harmonic oscillations. Indeed, 
in the particular case of strong collisions (Tl), Eq. ( 1181) reduces for 7 = j t to 

D = -4(<5 2 + ft 2 ) 3 < 0. (20) 

In the particular case of intense external field (T3), Eq. (fl~8l) reduces for ft ^> 74, 7 to 

J D = -fi 4 [4fi 2 -( 7 -7*) 2 ] <0. (21) 

Provided that 7 7^ j t , Table 1 and ensuing discussion at the end of section [2] suggest to 
measure 5 2 and ft 2 in the units of (7 — 7 4 ) 2 . Upon introducing a > and /3 > through 

5 2 = a(7-7*) 2 , fi 2 = /3( 7 -7*) 2 , (22) 

Eq. ([TBI becomes 



£ = -4(7 - 7t ) 



/3 3 + /3 2 ( 3a - - ) + /3a(3a - 5) + a{a + l) 2 



4, 

= -4( 7 - 7t ) 6 /i(a,0). (23) 

In the particular case of a zero detuning, 5 = a = (Torrey's case T2), the discriminant 
is greater than or equal to zero and all roots are real provided that < 1/4, or 

ft < b^pl . ( 24 ) 
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Obviously, each of the regions D = can be attained by a suitable choice of physical 
parameters (see figures [T] and [2]) in Torrey's case T2. Note in passing that h(a,(3) con- 
sidered as the cubic polynomial in /3 has only positive coefficients for a > 5/3. Thus any 
real root of h has to be negative. Consequently D(a, (3) < for any f3 > and a > 5/3. 
Actually a much stronger statement can be proven: 

D(a,P)< for any (3 > and a > 1/27. (25) 

The proof proceeds as follows. Considered as a cubic polynomial in 0, one has h(a, 0) > 
for any a > 0. Thus, as a straightforward consequence of the intermediate value theorem, 
h(a,(3) possesses a negative real root for any a > 0. Now it suffices to show that the 
discriminant dh of h(a,(3) considered as a cubic polynomial in is negative. This is 
indeed the case. Upon tedious but straightforward calculations one finds 

a 1 ix3 f < 0, a > 1/27, , , 

dh = ol (27a — 1 < _ _ I lnn 26 

h 16 v ; \ > 0, < a < 1/27. v ; 

The first inequality prohibits any additional real root of h(a,(3), and hence any root for 
f3 G (0, oo) and a > 1/27. Because h(a,0) > and h(a,f3 — > oo) — > +oo, one has 
necessarily h(a,f3) > for (3 £ (0, oo). The absence of zeros of D(a,(3) for a > 1/27 
is demonstrated in figures [T] and IU which plot the functional dependence of D c (a, (3) = 
—D/108 [see 122])] on (3 for selected values of a. Figure [2] shows the boundary of the D c < 
region in the (Q 2 ,5 2 ) plane. The boundary is exactly described by [28] 

a = -~ (2 + 3/3) + jj v 7 ! + 27/3 cos(# + 0„), (27) 



where 



= 1 cos" 1 



8 - 27/3(20 + 27£) 



(28) 



8(1 + 27/3) 3 / 2 

and = for the part of the boundary between the origin and the cusp point, whereas 
8o = 4:7i I '3 for the part of the boundary between the cusp point and (3 = 1/4. As it 
will be demonstrated in section [6] the critical value a = 1/27, which corresponds to the 
cusp point in figure [21 is related to the sufficient and necessary conditions ( 1501 for the 
occurence of a triply degenerate real root. Since the conditions for a triply degenerate real 
root can also be satisfied (see section [6] below), all functional forms of solutions of the 
Bloch equations listed in Eqs. (fTT|) are physically achievable. 



5 Decay constants 

In this section, various parameter ranges corresponding to each of the solutions in f }T7]) 
of the Bloch equations (T5]) are identified and decay constants of various solution types 
are explicitly determined. A prerequisite for that is a relation between roots of the cubic 
polynomial A(p) [see f lTU]) ] and the physical parameters of the Bloch equations The 



9 



latter is provided by means of Cardano's formula j2H 125] • Cardano's roots z% of a cubic 
polynomial f llUp have conventionally been written as follows 



where 



Zl 
Z3 

ft 



-g o 2 + (ft + ft), 

-i a 2 - ~ (5+ + ft) + \iV3(S + - ft) 

~02-| (5+ + ft) - iiV3(s+ - ft) 



R+VD C , S^ = \JR-^D C , 

and, on substituting the values of the coefficients aj according to Eqs. ( 1TT|) -( 1T3"|) . 

9a2ai — 27ao 



2a\ 



54 

3ai — a\ 1 



(7 - It) 



\ (A 2 " 25 2 ) - 1 (7 - 7 t ) 2 



i? 



Note in passing that «x in (jTTj) can always be obtained from the first Cardano's root (129]) . 



3 (S 2 + " g (7 - 7t) 2 , 
D 1 
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Q 3 + R 2 =■ a 

^ 108 108 

[4( 7 - 7t ) 4 5 2 + 4(5 2 + ft 2 ) 3 - ( 7 - 7t ) 2 (fi 4 + 205 2 fi 2 - 85 4 )] 



(29) 
(30) 

(31) 
(32) 

(33) 



1 S/^V^D 3 o 



(34) 



provided that (i) a cubic root y/x of a real number x is chosen to be a real number with 
the same sign as x (Rl); (ii) the respective cubic roots v^r ± of complex conjugate 
numbers x ± iy remain complex conjugate numbers (R2) [26]. As an example, in the 
particular case of strong collisions (Tl), one has a 2 = 37, R = [see (13T!) ]. 



^4 



±\JD C = ±v = ±-= V5 2 + ft 2 

v3 



(35) 



[see ( 13T)j) combined with that D c = — D/108], and the Cardano's formula ( 121?]) reproduces 
the known roots [21 [8] 



Xi = -7, 



z 2 3 = - 7 ± W5 2 + fi 2 . 



For 7 7^ 7t one finds 



ft 



(7 - 7t) p (/3-2q)-1±^ v / M^9)- 



(36) 



(37) 
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The real roots could be explicitly determined according to formula [25] 

x l+1 = -~ a 2 + 2Re (e l2wl/3 S) . (38) 

The root formula (1351) is fully compensated: it does not matter either which of S+ and S- 
has been taken for S, or which of the cubic roots has initially been taken for S± in (130]) . 
The set of real cubic roots given by (1551) remains invariant under any of the above choices, 
with a particular choice affecting only an irrelevant root permutation within the root set. 

In the particular case of a zero detuning, 5 = 0, or exact resonance (Torrey's case T2), 
A(p) factorizes into the product [see ([6])] 

A(p) = (p + 7t )[(p + T*)(P + 7) + (39) 

Thereby one needs only to solve a quadratic equation to obtain the roots (Appendix of 
Ref. [2J; section 3.5 of Ref. 0) 

xi = -It, z 2 , 3 = ± v/( 7 -7t) 2 -4fi 2 . (40) 



6 Degenerate real roots 

In virtue of the definition ( I33p . a necessary and sufficient condition for D > (or equiva- 
lently D c < 0), i.e., for purely damped nonos dilating solutions, is 

Q 3 < -i? 2 < 0, (41) 

which in turn requires Q < 0. Given ( |32|) . the latter implies 

5 2 + fi 2 <i( 74 - 7 ) 2 ^« + /3<i (42) 

which is necessary (but not sufficient) condition for D > 0. Obviously, the condition (142 p 
cannot be satisfied in a nonzero driving field for 7 = j t . 
One has D = [or D c = - see ([33])] if 

g 3 = -R 2 < 0. (43) 

Then according to (1301) 

5 ± = v^R, p c = 0). (44) 

Provided that the second inequality in ( 1431) is a sharp inequality, one has a doubly de- 
generate real root. Indeed, according to formula ( 135]) . A(p) has exactly two degenerate 
roots £2,3 = ^2,3 if and only if D = and R 7^ 0, and consequently, in virtue of the 
definition ( 133]) . Q < 0. A necessary condition for two degenerate roots is again given by 
f j4"2"]) . As an illustration, in the special case of 5 = (T2) the case of D = corresponds to 
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^ — It* — 7l/2 [see (|2~41) ] and R = (7 — 7i) 3 /216, which, according to our root convention 
Rl and R2, yields S± = = (7 — 7t)/6 [see 



^2,3 = g - , ^ 

and hence a doubly degenerate real root [see ( l40i) ]. 

The conditions for a triply degenerate real root can be satisfied for a nonzero Q only 
with a nonzero detuning. Indeed, according to (|38|) . a cubic polynomial can have a triply 
degenerate real root if and only if S± = 0, in which case the rotating term in (l38l) vanishes, 
and 

1 1 

21,2,3 = -7^2 = -3(7 + 27*) < 0. (46) 

According to (!30|) . S± = is possible if, in addition to D c = 0, one has also R = 0. The 
above conditions can only be satisfied if R = Q = [cf the first equality in (|33|) ]. or, given 
Eqs. (I3TI) and (1321) . when simultaneously 

2( 7 - lt f = 9fi 2 -185 2 , (47) 

(7~7*) 2 = 3(5 2 + fi 2 ). (48) 

In term of 6 and Q 

S 2 = ^(l-lt) 2 ^a = ± (49) 

fi 2 = 85 2 = |(7-7*) 2 ^/3 = ^, (50) 

and one finds a + j5 = 1/3 [see (142!) ]. Thus unless 7 = 7* (Tl) one can always attain the 
case of a triply degenerate real root for Q ^ 0. 

According to (|26l and figures [1] and [2] one has D c (a, (3) > for /3 > and a > 1/27. 
Therefore, the case of degenerate real roots is confined to rather small detunings \S\ < 
0, 19245(7— 7 t | [see (I2"2"j) ]. The dependence on (3 in figures [T] and [2] is extended to unphysical 
values of /3 < 0, which correspond to an imaginary magnitude of Q, in order to make the 
cubic dependence on (3 [see the square bracket in ff23l) ] transparent. For a 6 (0, 1/27), the 
zeros of D c (a,f3), considered as a cubic function of (3 > 0, occur in pairs (see also figure 
E]). Indeed, the coefficients aj of a general cubic polynomial such as that in ( TTUj) can be 
alternatively expressed in terms of cubic roots z\, z 2 , and z 3 as a x = ziz 2 + ziz 3 + z 2 z 3 and 
ao = —212:2-23. Given that the constant term of h satisfies ao > for a > 0, and combined 
with the existence of a negative real root of h, the real roots of h have to come necessarily 

either with the signs — h + or . The fact that in the case of h one has a\ < for 

a < 1/27 then excludes the option. 

According to figure El with increasing a: 

• (i) the interval between the pair of zeros of D c (a,(3) along a (3 trajectory decreases 
and 

• (ii) the interval middle point shifts slightly to larger values of [3 (see also figured]). 
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Since a < 1/27, and hence the condition (150]) is not satisfied, each of the zeros of the pair 
corresponds to a doubly degenerate real root. Along the boundary between the origin and 
the cusp point one has h = in fl37|) . the term (3 — 2a increases monotonically within the 
boundaries 

0</3-2«<^, (51) 

and the doubly-degenerate root changes monotonically between and — (7 + 27t)/3. 
Along the boundary between the cusp point and (3 = 1/4, the term 0— 2a in (137]) continues 
to increase monotonically within the interval 

— <f3-2a<-, (52) 

and the doubly-degenerate root changes monotonically between — (7 + 2j t )/3 and — (7 + 
7*)/2- 

The case of a triply degenerate real root, which according to ( 1501) occurs at a = 1/27 
[see (126]) ]. corresponds to the case when the interval between the pairs of doubly degenerate 
zeros of D c along a given /3-trajectory (parallel to the x-axis in figure |3]) reduces to zero. 
The latter corresponds to the cusp point of the boundary shown in figure EJ which separates 
the D c > and D c < regions in the (Q 2 , S 2 ) plane. 



7 Steady state solutions and remaining numerical con- 
stants 

Cubic roots of A(p) determine the decay rates in ( fl71) . The remaining part is to determine 
the numerical constants in ffTT]) . Two of the numerical constants could be determined from 
the initial conditions for the Bloch variables [see (IE 



A + A 1 + A 2 , D < 

. A + A x + A 2 + A 3 , D>0 

" l({, > \ A + A ± + A 2 , D = (double root) 

Aq + Ax, D = (triple root) 

and for the first derivative of the Bloch variables 

-KxAx - bA 2 + A 3 , D < 

, , -KxAx - k 2 A 2 - k 3 A 3 , D > ( . 

" l( ' ^ -k x Ax - k 2 A 2 + A 3 , D = (double root) 1 > 

-KxAx + A 2 , D = (triple root). 

In the latter case the left-hand side is provided by ([2]) taken at t — 0. A great deal of 
simplification can be achieved when some of the constants (e.g. a steady state solution Aq 
and Ax) could be determined in advance, before one makes use of the initial conditions. 
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Since A(p) has only nonzero roots, the steady state solution A can be determined from 
the PFD's listed by Eqs. ^ as 



A ° = 7777 lim n Pf(p) = — x { l^ltw eq , I = 2 (55) 
A(0) a ° ( 7(7 t 2 + I = 3 

where ao is given by ( 1T3|) and we have employed (J2J) in determining \im p ^ pf(p). Provided 
that the root x x = —k x is nondegenerate (e.g. D^O), the constant A x can be determined 
from the PFD's listed by Eqs. fTT5|) as [2] 

A 1 =^ b-*i)f(p) = f_M (56) 

p->*i A{p) (Xi - Z2)(Xi - z 3 ) 

Obviously, in the case D > all the constants Aj, j = 1, 2, 3 can be determined by a cyclic 
permutation of ( 156|) . For nearly degenerate roots one could instead use the expressions 

A 2 = — [x 3 (u -Ao- Ax) - (u' Q - XxAx)} , (57) 

x 3 - x 2 

A 3 = — [-x 2 (u - A - Ax) + (u' - x x A x )\ . (58) 

x 3 - x 2 

In the case of a triply degenerate real root, one determines A\ and A 2 straightforwardly 
from the initial conditions ( 153"]) and fl54l) and the knowledge of A [see fl55l) ]. A 3 could be 
determined from 

A 3 = hm — = f{-Kx). 59 

p-»~«i 

One finds that the substitution of £1,2,3 = — «i from (T4"6"j) . for p in ([2]) amounts to replacing 

2 1 

(p + 7) 3 (l-lt), (P + Jt) -»• -g (7-Tt)- ( 60 ) 

Thus in any case it is possible to determine v4 and one of the constants Ax, A 3 directly 
from the knowledge of the roots and x x = — «i of the product pA(p). 

So far we have ignored a nearly trivial case of 74 = 5 = 0, in which case ao = [see 
(fl3"|) ]. In general, one of the roots of A(p) is Ki = if a = [see (llOl)]. Then A(p) 
in ( jTUi) factorizes into a product of p and a quadratic polynomial, and all the roots of 
A(p) can be straightforwardly obtained. A necessary modification of the PFD, and of the 
recurrences for the coefficients Aj, to the case of a doubly degenerate real root p = is 
rather straightforward. Eqs. (l6Tj) are modified to 

{Ao 1 A 1 1 A 2 (p+b)+A 3 n n 

f+f+Jfe+F& ^>0 (61) 
f + # + + D = ( double root ) 
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Numerical constants are determined as follows. Instead of A by Eq. (|55p . one determines 

1 1 ( {-j6Qw eq ) , 1 = 1 

Al = XTnT lim n P 2 /(p) = x <^ 7 fi TtWeg , I = 2 (62) 

A(0)P "° Z2Z3 { lh? + 5 2 }w eq , 1 = 3 

Since it is no longer possible to obtain A through ( 15"5"|) . as the second constant for D > 
in (!6T|) one could determine 

^i(m + m = to(tm, (63) 



2 \ z z 

where z is one of the complex conjugate roots [discussed earlier in connection with (TIE]) ]. 
Eq. (|56|) applies for Aj only for j = 2, 3. In the case of a doubly degenerate real root 
x 2 = x 3 , 

A 3 = lim (P- ^) 2 /(P) = Z(£g). ( 6 4) 

We have excluded here the unrealistic undamped case of 7 = 74 = 0, which reduces to 
the precession of a classical gyromagnetic moment in a magnetic field [31 [8j , and which 
comprises the case a 2 = oq = leading to a triply degenerate root of A(p). 

In the absence of pwe phase relaxation processes (e.g. atoms in a dilute vapor cell; 
single molecule in solid hosts at superfluid helium temperatures [12]; localised surface 
plasmons), a pure dephasing rate is absent. Then the above expressions simplify according 
to the substitutions 7 —> 2j t and (7 — r y t ) 2 It- 



8 Discussion 

8.1 PFD over real number field 

Our treatment is unique in that it takes full advantage of the PFD over real number field, 
which made it possible to find and classify all analytic solutions. Surprisingly enough, 
the PFD has been nowhere mentioned in the context of Torrey's solution of the Bloch 
equations for a two-level system [21 EE El El ^9 anc ^ nas no ^ been used in its full generality. 
Surprisingly enough, earlier works [21 H] can be characterized by taking into account only 
one of the possible PFD's listed in Eqs. flm (see Eq. (41) of Ref. 0). 

Laplace transform combined with the PFD over complex number field has been em- 
ployed to solve related problems of the Bloch equations for a three-level system by Bernard 
et al |30j [see Eq. (4) therein] and five-level kinetics by De Vries and Wiersma [cf Eq. (A6) 
in Appendix A of Ref. [31]]. However, by making use of the PFD in the complex domain 
one can no longer guarantee that all the numerical constants in (16 7p are real numbers. The 
complex PFD obscures the fact that although any linear combination U = Xi^~ Klt +X2^~ K2t 
with real K\, k% > and Xu X2 7^ can be recast as V = £ie~ 6 * cosh(st) + ^2^ ht sinh(st) 
with some real b and s, the reverse is only possible for b ± s > 0. However, if the lat- 
ter condition holds, then V becomes an awkward recasting of U which obscures that V 
reduces to a sum of two simple exponential decays (cf Ref. [28]). 
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8.2 Earlier work 



Torrey [2] did not additionally make use of Cardano's formula. His general solution was 
confined to the case D < 0, when the characteristic determinant A(p) has a pair of complex 
conjugate roots. Therefore Torrey [2] missed the solutions corresponding to the parameter 
range D > 0. Hore and McLauchlan formally determined the cubic roots by Cardano's 
formula [4]. However, solution given by their Eq. (3) was limited to D > 0. Indeed, 
coefficients Aj given by their Eq. (4) become singular for degenerate roots. Furthermore, 
although the denominator in the expression for Ai in their Eq. (4) is formally correct for 
D > [cf our f l56|) ]. its numerator does not appear to be equal to the value of fi(—Ki) 
defined by Eqs. (Q. Additionally, Hore and McLauchlan solution misses a constant term 
[see Eq. (3) in Ref. [I] with the second functional form of solutions of the Bloch equations 
in our (fT7|) ]. 

Noh and Jhe [28] have correctly pointed at the incompleteness of Torrey's solution [2] . 
They employed a slightly asymmetric form of the Bloch equations, which resulted when 
u and v in ([3]) were defined without the pref actor of two. However, they did not solve the 
Bloch equations ab-initio and did not make use of the Laplace transform. Instead they 
employed a trial Ansatz and implicitly employed Cardano's formula to classify different 
solutions. Their solution is only limited to the special (although the most important and 
the most studied) case of wo = w eq = — 1. In their approach, Noh and Jhe [28] were 
required to look at the initial conditions for the second derivative of the Bloch variables, 
in spite that the Bloch equations ([2]) do only contain the first derivatives. Neither Noh 
and Jhe [28J made a connection between D c in Cardano's formula and the discriminant 
D of a cubic polynomial [see ( 13"3"|) and (JEJ)]. In the case of D > 0, Noh and Jhe solution 
is expressed in terms of hyperbolic sine and cosine (see Eq. (12) in Ref. [28]) - they 
overlooked that it can be simplified to a sum of three exponentially decreasing terms. 
Additionally, they did not identify different parameter ranges of D = in terms of the 
intrinsic parameters of the Bloch equations [see our (fl8|) . ( 142|) . (j4"4"|) ]. 

8.3 Transients and unorthodox solutions 

It has been shown that unorthodox solutions of the Bloch equations are confined to rather 
small detunings 5 2 < (7 — 7 t ) 2 /27 and small field strengths Q 2 < 8(7 — 7 t ) 2 /27, which 
are readily accessible experimentally. Figure H] shows that, regarding time dependence of 
the Bloch variables, there is a very smooth transition between the respective regions of 
D c < and D c > 0. Even if a r were an order of magnitude larger than the boundary 
value of a r ~ 0.35 for f3 = 0.2, the resulting time- dependence of w(t) would still resemble 
a featureless exponentially damped curve. The only indication that one is outside the 
D c < region is that the steady-state value of w becomes marginally smaller than the 
curve maximum. Therefore, in the parameter range of unorthodox solutions and within a 
substantially larger parameter subrange of Torrey's solutions proximal to the unorthodox 
solutions, the Bloch variables approach their respective steady states without reaching 
any significantly higher values in a transient region. 
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9 Conclusions 



A complete classification of unorthodox solutions of the Bloch equations for a two-level 
system in time-harmonic driving classical field, which do not have the form of familiar 
Torrey's exponentially damped harmonic oscillations, has been provided. Parameter range 
of the unorthodox solutions has been shown to be readily accessible experimentally. Time 
dependence of unorthodox solutions is characterized by rather featureless exponentially 
damped behaviour. The unorthodox solutions are essential for a reliable description of 
many different magnetic resonance [H El El HI El CEO EES] and quantum optics two-level 
systems (21 El El El QUI Ell EE21 EH EH]- The complete set of solutions could also provide 
a testing ground for general operator techniques involving exponential solutions of differ- 
ential equations for a linear operator [321 E21 EH ES] . A F77 code used to generate plots 
here is freely available [36J. 
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A Partial fraction decomposition over the field of real 
numbers 

Suppose there exist real polynomials f(x) and g(x) ^ 0, such that 

h(x) = ^ ' (65) 

By removing the leading coefficient of g(x), we may assume without loss of generality 
that g{x) is a polynomial whose leading coefficient is one (i.e. monic polynomial). By the 
fundamental theorem of algebra, we can write 

g(x) = {x- xi) jl ■■■(x- x m ) jm (x 2 + bix + Ci) kl ■ ■ ■ (x 2 + b n x + c n ) kn , (66) 

where Xi, . . . , x m , bi, . . . , b n , ci, . . . , c„ are all real numbers with b 2 —4ci < 0, and ji, . . . ,j m , 
ki, . . . , k n are positive integers. The XjS correspond to real roots of g(x), and the terms 
(x 2 + biX + [the so-called irreducible quadratic factors of g(x)] correspond to pairs of 
complex conjugate roots of g(x). The partial fraction decomposition of h(x) is 

f(x) m ii n ^ g x q. 

h{x) = W) = q{x) + ^ + ^ (x 2 +lx +ly ' (67) 

j l=l r=l \ %) l=l r= i \ I D 

Here q(x) is a (possibly zero) polynomial, and the A ir , B ir , Ci r , bi, and Cj are all real 
constants. A further information on the partial fraction decomposition can be found in 
Ref. [231. 



B Summary of inverse Laplace transform formulae 



C'^s- 1 } = 0(t), 



{{s + a) 2 + co 2 J 

cr x L * + ° 2 | = e(t)e-^cosM, 



1 1 „, , t 



n-l 



\(s + a) n ) w (n- 1)! 

£-i{ s - 2 } = Q(t)t, (68) 

where 0(t) is the Heaviside step function whose value is zero for negative argument and 
one for positive argument. 
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Figure captions 



Figure 1. - Plots of D c (a,f3) as a function of /3 for various values of the detuning S 2 
in the units of (7 — 7 t ) 2 /27. One has D c (a,j3) > for j3 > 0, provided that a r > 1 
(a > 1/27 « 0.037). 

Figure 2. - Zoom-out view of Figure 1 showing cubic curves of the plots of D c (a, (3) as 
a function of (3. 

Figure 3. - The boundary separating D c > and D c < regions in the (Q 2 , 5 2 )-plane. 
The cusp point corresponds to the triply-degenerate root. Along the boundary between 
the origin and the cusp point, the doubly-degenerate root changes monotonically between 
— 7t and — (7 + 27^/3, whereas between the cusp point and j3 = 1/4 the doubly-degenerate 
root changes monotonically between — (7 + 27 t )/3 and — (7 + 7t)/2, respectively. 

Figure 4. - Time dependence of of the Bloch variable w(t) for 7 = 0.4, 7 t = 0.1, and 
(5 = 0.2 for various values of the detuning 5 2 in the units of (7 — 7 t ) 2 /27. The first two 
smallest values of a r correspond to the D c < region (/3 = 0.2 corresponds the boundary 
value of a ~ 0.01299, or a r ~ 0.35) and the remaining two values of a r correspond to the 
Torrey (D c > 0) region. Note in passing a very smooth transition in the time dependence 
of w(t) between the respective regions of D c < and D c > 0. 
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